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R " radius  of  satellite 
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Ap  ” Debye  length 


“ satellite  potential 


Tg  - electron  temperature 


M+  **  Ionic  mass 


Vg  " satellite  velocity 


The  ratio  a is  a measure  of  the  amount  of  perturbation  created  by  the  motion 

of  the  satellite.  r P — ^ 

The  preliminary  results  obtained  jdemonstrate  in  a quantitative  manner 
the  dependence  of  a - f(altitude)  and  o - f(R^)  for:  24<RJ)<56  and  10<{*j^<18. 
We  find  that  a at  an  altitude  of  aboiit  300  km  is  about  2.5  to  3 orders  of 
magnitude  smaller  than  a at  an  altitude  of  about  500-600  km  and  a at  this 
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about  1100  Jb.  We  also  find  a for  IL  $54-56  to  be  3.5  order  of  magnitude 

*(„>«.<  P 

smaller  than  o for  R^  $ 24-30.  The  investigation  already  performed  yielded 
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I.  OBJECTIVE  OF  STUDY 

The  overall  objective  of  the  present  study  is  to  determine  the  angular 
distribution  of  the  ion  current  around  an  ionospheric  satellite  and  investigate 
phenomenologically  the  dependence  of  the  normalized  current  on  some  plasma 
parameters. 

The  overall  technological  - application  objective  is  to  determine  the 
extension  of  the  spatial  disturbances  created  by  the  motion  of  the  satellite 
and  their  interaction  with  thermal  measurements  performed  by  plasma  probes 
(e.g.,  current  collectors  of  various  kinds)  onboard  satellites.  From  the 
application  point  of  view  the  above  study  should  aid  in  assessing  the  validity, 
reliability  and  quality  of  particle  and  field  in-situ  measurements  and  help  in 
the  planning  of  Spacelab  experiments. 

Scientifically  the  study  should  yield  clues  as  to  the  existence  or 
nonexistence  plasma  oscillations  and  shock  patterns  ahead  and  behind  the 
satellite  for  specific  combinations  of  plasma  parameters. 

Technically,  and  from  the  application  point  of  view,  the  study  should 
have  a strong  bearing  on  the  planning  and  design  of  probe  packages  to  be 
mounted  on  future  satellites,  and  in  particular  on  the  Shuttle/Spacelab.  It 
should  be  realized  that  now  with  the  advent  of  Shuttle/Spacelab  many  ionospheric 
magnetospheric  and  spaceplasma  physics  experiments  are  planned.  These  types 
of  experiments  will  utilize  ensembles  of  diagnostic  probes  mounted  on  booms 
ejectable  instrument  packages  (known  as  throw  away  detectors) and  on  maneuverable 
subsatellites  of  varying  degrees  of  sophistication.  The  types  of  instruments 
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considered  at  present  for  such  experiments  include  low  energy  current 
collectors  such  as  Langmuir  probes  and  retarding  potential  analyzers 
(of  various  kinds)  and  antennas.  In  the  past,  perturbations  created  by  the 
motion  of  the  satellite  were  never  taken  too  seriously  in  as  far  as  their 
potential  contribution  to  the  reliability  and  validity  of  low  energy 
measurements  is  concerned.  This  situation,  however,  cannot  continue  in 
the  coming  Spacelab  era.  Therefore,  the  main  objective  of  the  present  study 
and  the  suggested  study  of  the  coming  few  years  is  mainly  technological/ 
application  oriented,  without  diminishing  the  value  of  the  scientific  outcome. 

In  the  last  contract  period  we  have  focused  mainly  on  the  investigation 
of  the  ratio  of  ion  current,  in  the  wake  as  a function  of  some  plasma  para- 
meters. In  other  words,  we  focused  on  investigating  the  variations  of 
the  ratio  of  the  current  in  the  wake  normalized  with  respect  to  the  current 
as  measured  when  the  probe  looked  ahead  (ma)  as  a function  of  various  plasma 
parameters,  and  in  particular,  the  ratio  of  the  radius  of  the  satellite 
divided  by  the  Debye  length. 
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II.  GENERAL  SCIENTIFIC  AND  TECHNOLOGICAL  BACKGROUND 
- BRIEF  DISCUSSION 


It  is  well  known  that  a spacecraft  orbiting  around  the  earth  or  orbiting 
in  any  other  planetary  environment  interacts  with  its  space  environment  thereby 
causing  significant  perturbations  in  the  near  and  far  vicinity  to  its  surface. 
The  interaction  between  the  spacecraft  and  its  environmental  plasma  is  mutual 
in  that  the  phenomena  involved  are  coupled  by  the  effects  on  the  space  plasma 
and  on  the  spacecraft. 

The  spacecraft  itself  acquires  a potential  due  to  the  collection  of  the 
surrounding,  charged  particles,  and  due  to  other  charging  mechanisms,  such  as: 
photo-emission,  secondary  particle  emission,  magnetic  field  effects,  energetic 
particle  bombardment,  etc.  (e.g, , Whipple,  1965;  Samir  and  Willmore,  1965; 
Whipple,  1977;  Samir  and  Willmore,  1966).  The  interaction  implies  that  there 
is  a transfer  of  charge  from  the  environment  to  the  spacecraft  and  from  the 
spacecraft  to  the  environment.  In  stating  the  above  we  ignore  the  transfer 
of  charge  between  different  elements  on  the  spacecraft  surface;  namely,  we 
assumed  the  surface  to  be  a perfect  conductor.  Needless  to  say,  that  in 
practice  this  is  usually  not  the  case. 

The  surrounding  space  plasma  is  strongly  perturbed  by  the  satellite 
motion  and  a wake  zone  depleted  unequally  of  ions  and  electrons  is  therefore 
formed  behind  the  spacecraft.  In  addition  to  the  creation  of  the  wake,  the 
spacecraft  excites  plasma  oscillations  in  the  electrostatic  ion  plasma  mode 
(Samir  and  Willmore,  1965).  A potential  well  behind  the  spacecraft  at  some 
distance  from  its  surface  downstream  should  also  exist.  It  can  be  antici- 
pated (Liu,  1969;  Gurevich  £t  aJ. , 1970;  Gurevich  and  Pitaevski,  1975;  and 

) 


3 


03  02  019 


T 


I 


Liu,  1975)  Chat  plasma  oscillations  and/or  particle  energization  mechanisms 
are  associated  with  the  potential  well.  Moreover,  the  fact  that  there  is  a 
steep  gratient  in  electron  and  ion  densities  between  the  wake  zone  and  else- 
where around  the  satellite  gives  rise  to  the  prediction  that  some  kind  of 
plasma  instability  should  exist  in  the  disturbed  zones. 

Mathematically,  the  complex  of  phenomena  involved  in  the  interaction 
between  a spacecraft  and  its  surrounding  space  plasma  has  to  be  described  in 
a self-consistant  manner.  This  means  that  the  Vlasov  equation  written  for 
each  charged  component  of  the  plasma  together  with  the  Poisson  equation  have 
to  be  solved  simultaneously  in  a self-consistent  way  subject  to  realistic 
boundary  conditions.  It  is  well  known  that  this  is  a difficult  task  and  no 
general  solutions  yet  exist  despite  the  fact  that  various  aspects  of  the 
interaction  were  studied  theoretically  for  over  a decade.  For  general  review 
papers  on  the  subject  we  sight  Liu,  1969;  Gurevich  e£  al . , 1970;  Gurevich  and 
Pitaevski,  1975;  Samir,  1973;  Whipple,  1977.  The  available  wake  models  require 
confirmation  or  varification  via  comparing  theory  with  experiment.  At  the 
present  time  the  status  of  such  verification  or  confirmation  is  not  satis- 
factory. The  few  attempts  made  to  compare  results  from  in  situ  measurements 
with  theoretical  sheath  and  wake  models  are  to  be  considered  preliminary  at 
best.  This  situation  is  due  largely  to  the  small  amount  of  relevant  in  situ 
information  available  for  such  an  analysis.  It  should  also  be  realized  that 
the  in  situ  information  available  is  often  fragmentary  in  that  not  all  the 
physical  parameters  required  to  carry  out  a meaningful  theory  experiment 
comparison  are  available.  The  present  status  of  the  study  is  therefore 
exploratory  and  there  is  a great  need  to  examine  larger  samples  of  electron 
and  Ion  densities  (currents)  and  potential  field  measurements  relevant  to  the 

i 

study  of  spacecraft  space  plasma  interaction. 
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There  can  be  little  doubt  that  perturbations  created  by  the  interaction 
have  profound  relevance  to  the  realiability  and  quality  of  particle  and  field 
measurements  made  from  satellites  and  rockets. 

The  study  of  available  in  situ  measurements  is  expected  to  help  in  the 
development  and  planning  of  future  active  and  passive  experiments.  It  should 
be  realized  that  the  planning  of  experiments  relevant  to  spacecraft  space 
plasma  interactions  is  now  in  progress.  The  planning  aims  at  utilizing  the 
shuttle  spacelab  facility  as  a near  earth  plasma  laboratory.  Until  now  we 
focussed  mainly  on  the  application  aspect  of  our  study.  However,  as  mentioned 
earlier  we  would  like  to  utilize  the  results  so  as  to  obtain  scientific 
information  relevant  to  future  studies.  Namely  we  would  like  to  perform 
a theory  experiment  comparison  via  comparing  results  from  the  measured 
ion  depletions  in  the  wake  with  (as  far  as  possible)  the  theoretical  computa- 
tions for  the  similar  plasma  parameters.  Hence,  examine  the  real  insignificance 
of  the  physical  assumptions  used  in  constructing  the  theoretical  models. 

A brief  general  discussion  regarding  theoretical  studies  is  at  this  point 
appropriate.  As  mentioned  earlier  the  problem  of  the  theoretically  calculated 
structure  of  the  disturbed  plasma  around  a moving  body  in  space  involves  the 
solution  of  a complicated  system  of  coupled  non-linear  partial  differential- 
integral  equations.  The  equations  consist  of  Vlasov  equations  for  the  ions 
and  electrons  and  the  Poisson  equation  relating  the  electric  field  to  the 
distributions  of  ions  and  electrons.  In  cases  of  stationary  bodies  as  well 
as  in  cases  of  moving  bodies  combinations  of  numerical  techniques  (e.g.  finite 
differences,  iterations,  etc.)  are  required. 

Various  approaches  have  been  used  in  the  past  to  deal  with  this  problem; 
however,  due  to  the  difficulties  involved  in  the  self-consistent  solution, 
simplified  assumptions  had  to  be  used.  Among  the  customary  ones  we  cite: 
ignoring  the  geomagnetic  field  influence  on  particle  collection,  using  simplified 


geometries  for  body  configuration,  using  simplified  surface  reactions,  assuming 
that  the  satellite  is  a conducting  body, and  assuming  that  it  is  justified  to 
treat  the  problem  as  a steady  state  one.  One  of  the  more  io.cent  theoretical 
studies  which  adopted  a self-consistent  mathematical  treatment  is  that  of 
Parker  as  j,iven  in  Parker,  1976.  Parker  adopted  some  of  the  simplifying  assump- 
tions but  has  advantages  over  other  theoretical  treatments.  Parker's  method 
uses  a numerical  grid  (or  mesh)  of  discrete  points  in  space  where  the  potential 
and  density  distributions  are  defined.  The  Poisson  and  Vlasov  equations  are 
represented  in  finite  difference  form  at  these  grid  points.  In  his  most  recent 
work,  Parker  used  a pillbox  shaped  body  where  the  surface  of  the  pillbox  coincides 
with  certain  rows  and  columns  of  grid  points.  One  advantage  is  that  the  grid 
points  are  unequally  spaced  so  that  a higher  density  of  points  can  be  used  near 
the  surface  and  a lower  density  of  points  used  further  away  from  the  surface. 

This  allows  a given  number  of  grid  points  to  be  used  efficiently.  Despite 
the  fact  that  this  model  seems  to  us  to  be  the  most  sophisticated  and 
realistic  wake  model  fundamental  questions  still  exist.  For  example,  the 
question  of  whether  or  not  is  it  justified  to  use  a kinetic  approach  at  certain 
distances  surrounding  the  body  and  a fluid-like  approach  at  distances  further 
away  from  the  body.  Or  is  it  justified  to  use  a fluid-like  approach  for 
specific  plasma  and  body  parameters  and  a kinetic  approach  for  other  sets  of 
parameters.  This  is  undoubtedly  a very  interesting  scientific  problem  which 
should  have  far  reaching  applications  in  cosmic  electro  dynamics  when  interest 
is  focussed  on  the  interaction  between  the  solar  wind  around  planets  and  the 
natural  satellite  and  interstellar  plasma  interacting  with  bodies  in  deep 
space.  For  details  on  the  difference  between  the  approach  adopted  by  Parker 
as  compared  with  approaches  adopted  by  other  authors  reference  is  given  to 
Parker,  1976. 
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In  the  present  study  we  were  concerned  with  one  aspect  only,  namely  with 


examining  the  amounts  of  ion  depletion  in  the  wake  of  the  S3-2  satellite  as 
a function  of  several  plasma  parameters,  since  we  believed  that  the  amount 
of  ion  depletion  in  the  wake  serves  as  a measure  of  the  significance  of  the 
disturbances  created  by  the  motion  of  the  satellite  to  its  environment.  We 
will  also  (at  a later  stage)  attempt  a comparison  between  the  measured  amounts 
of  depletion  with  results  obtained  from  Parker's  and  other  theoretical  models 
in  order  to  achieve  the  more  scientific  goal  of  assessing  the  relative  impor- 
tance of  assumptions  used  in  various  wake  models. 

In  the  next  section  we  present  and  discuss  some  preliminary  results. 
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III.  PRELIMINARY  RESULTS  - A DISCUSSION 


1.  Introduction 


During  the  last  contract  period  we  have  focussed  on  getting  the  data 
in  a form  that  lends  itself  best  to  our  purposes.  We  should  state  that 
we  have  encountered  some  unexpected  difficulties  in  the  data-handling  but 
we  do  have  now  samples  of  data  that  were  analyzed  yielding  some  interesting 
perliminary  results. 

The  main  research  effort  focused  on  extracting  selected  samples  of 
data  from  the  S3-2  satellite  in  order  to  obtain: 


(i)  the  variation  of  I+  = f(0),  where:  I+  = ion  current  and 
0 = angle  between  the  normal  (to  the  probe)  and  the  sat- 
ellite velocity  vector 

(ii)  the  basic  plasma  properties;  namely,  electron  temperature 
(T  ),  electron  density  (N  ),  plasma  potential  (<|>s ) (or: 
potential  of  the  sphericaf  electron  probe  with  respect  to 
the  plasma),  required  in  order  to  compute  the  following 
plasma  parameters: 


Note  that  M+  is  assessed  from  ionospheric  models. 
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In  addition  we  devoted  a considerable  effort  to  modifying  the  Parker 
wake/sheath  program  (e.g.  Parker,  1976)  trying  to  make  it  applicable  to 
the  case  of  a satellite  moving  in  the  terrestrial  ionosphere-magnetosphere 
plasma.  The  objective  of  the  latter  effort  is  (eventually)  to  be  able  to 
perform  a theory-experiment  comparison  which  is  after  all  the  best  test 
for  the  validity  of  theoretical  models  and  for  the  validity  of  the  physical 
assumption  used  in  the  construction  of  the  models. 


About  the  Data 

At  the  present  stage  of  the  study  we  have  used  ion-current  data  from 
several  selected  passes/Orbit  1131,  orbit  data  2/24/76  using  sensors  1-8. 

As  will  be  seen  in  the  next  section  the  quality  of  the  data  differs 
for  different  sensors. 

At  the  present  time  we  have  plots  showing:  I+  (in  amps)  = f (G.M.T.) 
and  we  have  the  angles:  9,  (n ,B) , (Vs,ji)  for  the  corresponding  times. 
Therefore,  we  are  able  to  examine  a reasonable  sample  of  I = f(0)  in  the 
altitude  range:  ^ 1100  km  to  - 300  km  (or  less).  The  quality  of 

the  data  varies  and  we  encounter  difficulties  in  the  perigree  part  of  the 
orbit.  This  is  unfortunate  since  the  information  at  the  low-altitude 
range  is  of  greater  scientific  value  to  us.  On  the  other  hand  even  the 
above  data  yielded  upper  limits  to  the  ratio: 


a 


I+  (wake) 
I+  (front) 


which  are  not  yet  known  or  published  in  the  literature.  Moreover,  from 
the  application  point  of  view,  the  latter  information  is  useful. 

It  should  be  realized  that  at  the  present  time  we  are  not  attempting 
to  computerize 


-r  ■+***  •'  **  ‘ fr  . 


in  a 'blind'  manner  but  rather  use  computer  outputs  together  with  a close 
'visual  inspection'  of  the  data.  Due  to  the  quality  of  the  data  this  seems 
to  us  a reasonable  procedure  to  follow.  This  approach  may  be  altered  in 
the  future.  It  is  our  experience  that  such  an  approach  'pays  off'  eventually 
although  it  is  slow  in  time.  It  should  be  possible  in  principle  and  if 
enough  resources  are  available  to  write  a pattern-recognition  program  to 
compute  [I+  (wake)]  but  we  doubt  the  usefulness  of  such  an  approach  when 
the  study  is  problem-oriented  using  relatively  small  samples  of  selected 
data.  After  our  preliminary  examination  of  the  selected  passes  from  orbit 
1131  is  completed,  we  will  use  more  measurements  from 
other  orbits  and  extend  the  study  to  utilization  of  E-field  data,  i.e. 

Michael  Smiddy's  data. 

In  addition  to  the  above  information,  we  have  computed  the  values  of 
the  parameters:  R^,  S and 
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2.  Presentation  of  Preliminary  Results 


Table  1 gives  the  ion  currents  (I+)  the  angles  of  attack  (0)  and  the 

ratios  (a)  for  the  measurements  from  sensor  1.  The  angle  6 . shows  the 

min 

angle  closest  to  0 + 0°  (which  is  the  'exact'  ram  direction)  and  I,  (0  . ) 

+ min 

shows  the  ion  current  as  measured  by  sensor  1 for  this  specific  angle.  To 

obtain  the  correct  I+  (0^^)  the  plots  mentioned  earlier  (i.e.  I+  = f (G.M.T. )) 

and  the  detailed  output  from  the  tapes  (in  table  form)  were  used.  The  same 

goes  for  0 and  I (0  ) except  that  for  the  maximum  wake  position  0 -*■ 

max  t max 

180°.  The  column  on  the  right  shows  the  corresponding  altitude.  The  column 
which  gives  Time  (G.M.T.)  is  for  our  purposes  of  comparing  information  from 
plots  with  information  from  tables,  and  is  of  no  value  to  the  reader. 

Table  2 also  for  sensor  1 provides  the  computed  values  of  the  nondimen- 
sional  parameters: 


, _ e*s 

" KT 


PN 


as  well  as  the  measured  values  of  the  plasma  properties: 


Ng  “ Density  (“electron  density) 

Tg  » Temp,  (-electron  temperature) 

PS  - <t>g  (=plasma  potential) 

VS  - Vg  (^satellite  velocity) 

together  with  the  corresponding  values  of  a and  H (^altitude). 
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TABLE  2: 

Plasma  Pa/ 

ameters 

Sensor 

J. 

Density 

PS 

Temp. 

VS 

RD 

PN 

a 

Ho  (km) 

0.377xl04 

-1.17 

1006. 

7.299 

27.86 

-13.52 

3.7xl0-2 

993.85 

0.367xl04 

-1.17 

1006. 

7.311 

26.25 

-13.52 

3. 2x10* 2 

981.28 

0.342xl04 

-1.34 

1180. 

7.324 

24.67 

-13.20 

3.2xl0~2 

968.69 

0.368xl04 

-1.22 

1011. 

7.337 

28.33 

-14.11 

3.1xl0"2 

956.06 

0.405xl04 

-1.19 

963. 

7.349 

29.46 

-14.37 

2.7xl0~2 

943.41 

0.411xl04 

-1.19 

963. 

7.362 

29.35 

-14.37 

2.3xl0-2 

930.72 

0.414xl04 

-1.19 

963. 

7.375 

30.05 

-14.37 

2.4xl0-2 

918.01 

0.437xl04 

-1.09 

976. 

7.387 

30.63 

-12.99 

2. 8x10* 2 

905.26 

0.427xl04 

-1.04 

1053. 

7.400 

28.39 

-11.50 

2.3xl0"2 

892.46 

0.437xl04 

-1.01 

1092. 

7.413 

28.72 

-10.75 

2.8xl0-2 

879.64 

0.450xl04 

-1.09 

1026. 

7.426 

30.85 

-12.35 

2. 6x10* 2 

866.81 

0.480xl04 

-1.09 

1026. 

7.439 

31.35 

-12.35 

2.4xl0“2 

853.99 

0.481xl04 

-1.09 

1026. 

7.452 

31.54 

-12.35 

2.2xl0~2 

841.16 

0.456xl04 

-1.08 

• 1032. 

7.465 

30.17 

-12.16 

2.1xl0~2 

828.34 

0.423xl04 

-1.07 

1036. 

7.478 

29.28 

-12.01 

2. 4x10* 2 

815.51 

0.410xl04 

-1.08 

954. 

7.491 

30.04 

-13.16 

2.4xl0~2 

802.65 

0.391xl04 

-1.08 

954. 

7.504 

29.34 

-13.16 

1.8x10* 2 

789.83 

0.392xl04 

-1.13 

963. 

7.517 

28.83 

-13.64 

2. 5x10* 2 

777.04 

0.393xl04 

-1.13 

963. 

7.531 

29.54 

-13.64 

2.0xl0~2 

764.26 

0.406xl04 

-1.13 

963. 

7.544 

29.61 

-13.64 

1.9xl0-2 

751.53 

0.432xl04 

-1.08 

1207. 

7.557 

27.42 

-10.40 

1.9x10* 2 

738.83 

0.472xl04 

-1.08 

1207. 

7.570 

28.66 

-10.40 

1.8xl0-2 

726.16 

0.551X104 

-1.19 

1123. 

7.583 

32.10 

-12.32 

1.7xl0-2 

713.54 

0.689X104 

-1.21 

1185. 

7.596 

34.94 

-11.89 

1.69x10* 2 

700.98 

0.822xl04 

-1.24 

1308. 

7.609 

36.50 

-11.02 

1.48xl0*2 

688.46 
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TABLE  2:  Plasma  Parameters  (cont)  Sensor  1 


Density 

PS 

Temp. 

VS 

KD 

PN 

a 

Ho (km) 

0.968xl04 

-1.24 

1308. 

7.622 

39.41 

-11.02 

1 . 33xl0~2 

676.01 

0.981xl04 

-1.24 

1308. 

7.635 

39.68 

-11.02 

1.16xl0-2 

663.65 

0.960xl04 

-1.32 

1254. 

7.648 

40.04 

-12.24 

1.16xl0'2 

651.38 

0.955xl04 

-1.35 

1302. 

7.660 

39.38 

-12.11  ' 

1.08xl0_2 

639.21 

0.954xl04 

-1.47 

1468. 

7.673 

36.94 

-11.64 

9.82xl0'3 

627.12 

0.981xl04 

-1.20 

1116. 

7.686 

43.83 

-12.57 

9.61xl0-3 

615.08 

0.949xl04 

-1.17 

1072. 

7.699 

43.30 

-12.69 

8.80xl0~3 

603.11 

0.959xl04 

-1.17 

1072. 

7.712 

43.9 

-12.69 

7.92x10' 3 

591.27 

0.916xl04 

-1.21 

1114. 

7.724 

42.36 

-12.66 

8.53xl0~3 

579.56 

0.930X104 

-1.55 

1449. 

7.736 

36.69 

-12.44 

8. 3xlO-3 

567.99 

0.967xl04 

-1.37 

1415. 

7.748 

38.02 

-11.26 

8. 32xlO_3 

556.54 

0.942xl04 

-1.32 

1405. 

7.761 

37.51 

-10.92 

8.89xl0~3 

545.21 

0.953xl04 

-1.48 

1278. 

7.773 

39.55 

-13.47 

9 -Oxio' 3 

533.96 

0.960xl04 

-1.48 

1278. 

7.785 

39.70 

-13.47 

8.78x10" 3 

522.87 

0.992xl04 

-1.48 

1278. 

7.797 

40.34 

-13.47 

8. 52xlo'3 

511.94 

l.OlxlO4 

-1.44 

1367. 

7.808 

39.82 

-12.28 

8.10x10" 3 

501. J7 

1.05xl04 

-1.42 

1412. 

7.820 

39.46 

-11.69 

7. 49xl0~3 

490 . 56 

l.lOxlO4 

-1.56 

1260. 

7.831 

42.70 

-14.40 

6.87xlo"3 

480.11 

1.20xl04 

-1.56 

1260. 

7.842 

44.67 

-14.40 

5. 86xl0"3 

469 . 79 

1.18xl04 

-1.87 

1637. 

7.853 

38.98 

-13.28 

-3 

4.76x10 

459.66 

1.13xl04 

-1.87 

1637. 

7.864 

38.14 

-13.28 

3. 26xl0~3 

449.71 

1.13xl04 

-1.87 

1637. 

7.875 

38.07 

-13.28 

2.73xl0'3 

439.95 

1.05xl04 

-1.70 

1223. 

7.885 

42.51 

-16.16 

2 .Oxio"3 

430.38 

1.12xl04 

-1.70 

1223. 

7.896 

43.87 

-16.16 

2 .OxlO~3 

421.00 

1.14xl04 

-1.99 

1535. 

7.906 

39.43 

-15.07 

1. 71xl0~3 

411.80 
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TABLE  2: 

Plasma  Parameters  (cont) 

Sensor 

_1 

Density 

PS 

Temp. 

VS 

RD 

PN 

a 

Ho  (km) 

1.16xl04 

-1.99 

1434. 

7.916 

41.52 

-16.27 

9.3xl0~4 

402.79 

1.17xl04 

-2.00 

1308. 

7.926 

43.38 

-17.78 

7.2xl0-4 

393.79 

1.27xl04 

-2.00 

1308. 

7.936 

45.11  ' 

-17.78 

4.2xl0-4 

385.42 

1.40xl04 

-2.00 

1308. 

7.945 

47.43 

-17.78 

3. 3xl0_4 

377.06 

1.51xl04 

-1.95 

1564. 

7.954 

44.98 

-14.50 

2.45xl0-4 

368.93 

1.56xl04 

-2.10 

1732. 

7.963 

43.75 

-14.16 

1.22xl0~4 

361.02 

1.56xl04 

-2.40 

2069. 

7.972 

39.78 

-13.49 

l.lSxlO-4 

353.31 

1.50xl04 

-1.92 

1609. 

7.981 

44.18 

-13.88 

345.84 

1.73xl04 

-1.92 

1609. 

7.989 

47.52 

-13.88 

1.26xl0-4 

338.60 

1.86xl04 

-1.92 

1609. 

7.997 

49.21 

-13.88 

8.4xl0~5 

331.56 

2.04xl04 

-1.98 

1791. 

8.004 

50.11 

-13.12 

4 . 76xl0~5 

324.78 

2.17xl04 

-2.19 

2430. 

8.012 

43.34 

-10.48 

318.25 

2.13xl04 

-1.87 

1690. 

8.019 

52.84 

-13.08 

312.00 

2.27xl04 

-1.84 

1597. 

8.026 

54.63 

-13.40 

5.29xl0-5 

306.00 

2.24xl04 

-1.91 

1566. 

8.033 

54.87 

-14.18 

1.8xl0-5 

300.25 

2.23X104 

-1.91 

1566. 

8.039 

54.69 

-14.18 

294.74 

2.25xl04 

-1.91 

1566. 

8.045 

54.87 

-14.18 

289.48 

1.89xl04 

-1.86 

1678. 

8.051 

48.52 

-12.93 

284.50 

1.95xl04 

-1.85 

1710. 

8.056 

48.91 

-12.58 

4.0xl0~5 

279.78 

1.90xl04 

-1.93 

1476. 

8.062 

51.90 

-15.20 

275.34 

1.84xl04 

-1.93 

1476. 

8.067 

51.20 

-15.20 

271.17 
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Table  3:  (cont'd) 


TABLE  4:  Plasma  Parameters 


Sensor  5 


Density 

PS 

Temp. 

VS 

RD 

PN 

I+(0max) 

l+(0min) 

Ho (km) 

4.11xl03 

-1.17 

1006. 

7.  306 

29.27 

-13.52 

3.56xl0~2 

987.24 

3. 86x10 3 

-1.34 

1180. 

7.318 

26.20 

-13.20 

3.16xl0~2 

974.66 

3.99xl03 

-1.34 

1180. 

7.  331 

26.63 

-13.20 

3.45xl0~2 

962.05 

4.46xl03 

-1.19 

963. 

7.343 

31.17 

-14.37 

2 . lxio”2 

949.41 

4. 70x10 3 

-1.19 

963. 

7.356 

32.01 

-14.37 

936.73 

4. 55x10 3 

-1.19 

963. 

7.  369 

31.49 

-14.37 

2.0xl0'2 

924.03 

4. 98x10 3 

-1.09 

976. 

7.381 

32.74 

-12.99 

2.4xl0~2 

911. 30 

4.99xl03 

-1.09 

976. 

7.394 

32.78 

-12.99 

2.4xl0”2 

898.54 

4. 56x10 3 

-1.01 

1092. 

7.407 

29.61 

-10.75 

2.3xl0-2 

885.71 

4.75xl03 

-1.04 

1070. 

7.420 

30.47 

-11.28 

2.53xl0~2 

872.89 

5 .03xl03 

-1.09 

1026. 

7.433 

32.09 

-12.35 

2.32xl0"2 

860 . 06 

5. 37x10 3 

-1.09 

1026. 

7.446 

33.16 

-12.35 

2.00xl0"2 

847.24 

5. 46x10 3 

-1.09 

1026. 

7.459 

33.42 

-12.35 

1.95xl0"2 

834.41 

4.88xl03 

-1.07 

1036. 

7.472 

31.45 

-12.01 

2.47xl0"2 

821.59 

4. 54x10 3 

-1.07 

1018. 

7.485 

30.55 

-12.27 

2.49xl0"2 

808.74 

4.46xl03 

-1.08 

954. 

7.498 

31.35 

-13.16 

2.46xl0~ 2 

795.90 

4.44xl03 

-1.12 

962. 

7.511 

31.13 

-13.59 

2.37xl0-2 

783.09 

4.27xl03 

-1.13 

963. 

7.524 

30.53 

-13.64 

2 . 31xl0~2 

770.31 

4.37xl03 

-1.13 

963. 

7.538 

30.86 

-13.64 

2 . 10xl0~2 

757.56 

4.57xl03 

-1.12 

990. 

7.551 

31.55 

-13.28 

1.91xl0-2 

744.84 

4.95xl03 

-1.08 

1207. 

7.564 

29.34 

-10.40 

1.92xl0-2 

732.16 

5.87xl03 

-1.16 

1141. 

7.577 

33.13 

-11.89 

1.64xl0“2 

719.51 

6.72X103 

-1.19 

1123. 

7.590 

35.44 

-12.32 

1.67xl0-2 

706.92 

8.25xl03 

-1.24 

1308. 

7.603 

36.38 

-11.02 

1.48xlO-2 

794.39 

9. 49x10 3 

-1.24 

1308. 

7.616 

39.72 

-11.02 

1. 19xl0~2 

681 . 89 
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TABLE  4:  Plasma  Parameters  (cont)  Sensor  5 


Density 

PS 

Temp. 

VS 

RD 

PN 

I+(0max) 

l+(0min) 

Ho (km) 

1.05xl04 

-1.24 

1308. 

7.629 

41.12 

-11.02 

l.OlxlO-2 

669.49 

1.14xl04 

-1.29 

1272. 

7.642 

42.07 

-11.83 

l.lxlO”2 

657.18 

1.04xl04 

-1.32 

1254. 

7.654 

41.78 

-12.24 

l.OxlO-2 

644.97- 

1.28xl04 

-1.47 

1468. 

7.667 

38.34 

-11.64 

l.OxlO-2 

632.83 

1.06xl04 

-1.47 

1468. 

7.680 

38.84 

-11.64 

9.5xlO~3 

620.78 

1.02xl04 

-1.17 

1072. 

7.693 

44.76 

-12.69 

8.3xl0-3 

608.76 

l.OSxlO4 

-1.17 

1072. 

7.705 

45.32 

-12.69 

7.27xl0~3 

596.86 

9.98xl03 

-1.17 

1072. 

7.718 

44.22 

-12.69 

6.9xl0"3 

585.09 

9. 7 3x10 3 

-1.55 

1449. 

7.730 

37.54 

-12.44 

7.3xl0-3 

573.46 

9. 64x10 3 

-1.55 

1449. 

7.743 

37.36 

-12.44 

7.1xl0-3 

561.95 

9. 4 2x10 3 

-1.32 

1405. 

7.755 

37.52 

-10.92 

6.6xl0-3 

550.56 

9.41xl03 

-1.  39 

1349. 

7.767 

38.38 

-12.05 

5.8xlO~3 

539.27 

9. 4 JxlO  3 

-1.48 

1278. 

7.779 

39.  36 

-13.47 

5.5x10' 3 

528.10 

9.57x10 

-1.48 

1278. 

7.791 

39.65 

-13.47 

517.09 

9.  78x10 3 

-1.48 

1278. 

7.803 

40.08 

-13.47 

5.6xlO~3 

506.25 

9. 80x10 3 

- 1.42 

1412. 

7.814 

38.18 

-11.69 

4.6xlO_3 

495.56 

1.05xl04 

-1.47 

1361. 

7.826 

40.24 

-12.59 

3 . 8xl0_  3 

485.04 

l.lOxlO4 

-1.56 

1260. 

7.837 

42.72 

-14.40 

3.7xlO-3 

474.66 

1.12X104 

-1.87 

1637. 

7.848 

37.96 

-13.28 

3.0xl0~3 

464.43 

1.14xl04 

-1.87 

1637. 

7.859 

38.15 

-13.28 

2.6xlO~3 

454.40 

l.lxlO4 

-1.87 

1637. 

7.870 

37.51 

-13.28 

2 . 36xl0~  3 

444.55 

1.12xl04 

-1.83 

1545. 

7.880 

38.75 

-13.92 

1 . 8xl0~  3 

434.89 

1.06X104 

-1.70 

1223. 

7.891 

42.64 

-16.16 

1. 3xlO~3 

425.42 

1.03X104 

-1.96 

1500. 

7.901 

37.60 

-15.19 

7.0xl0~4 

416.15 

1.13xl04 

-1.99 

1535. 

7.912 

?9.28 

-15.07 

5. 7xlO~4 

407.03 

l.lOxlO4 

-2.00 

1308. 

7.922 

41.93 

-17.78 

5.4xl0"4 

398.13 

24 
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TABLE  4:  Plasma  Parameters  (cont) 


Sensor  5 


Density 

PS 

Temp. 

VS 

1.21xl04 

-2.00 

1308. 

7.932 

1.31xl04 

-2.00 

1308. 

7.941 

1. 39xl04 

-1.96 

1507. 

7.950 

1.51xl04 

-1.95 

1564. 

7.959 

1.44xl04 

-2.40 

2069. 

7.968 

1.63xl04 

0 

CNJ 

1 

2069. 

7.977 

1.58xl04 

-1.92 

1609. 

7.985 

1.64xl04 

-1.92 

1609. 

7.993 

1.87xl04 

-1.92 

1609. 

8.001 

2.00xl04 

-2.10 

2430. 

8.008 

2.18xl04 

-2.19 

2430. 

8.016 

2.20xl04 

GO 

•H 

1 

1597. 

8.023 

2.17xl04 

-1.89 

1580. 

8.030 

2.22xl04 

-1.91 

1566. 

8.036 

2.22X104 

-1.91 

1566. 

8.042 

2.10X104 

-1.91 

1566. 

8.048 

1.68xl04 

-1.85 

1710. 

8.054 

1.81xl04 

-1.89 

1606. 

8.059 

lr81xl04 

-1.93 

1476. 

8.064 

1.35xl04 

-2.37 

2770. 

8.069 

RD 

PN 

l+(0max) 
1+ (0min) 

Ho (km) 

43.97 

-17.78 

389.45 

45.82 

-17.78 

380.99 

43.76 

-15.23 

372.75 

45.01 

-14.50 

364.74 

38.18 

-13.49 

356.93 

40.61 

-13.49 

349 . 35 

45.41 

-13.88 

342.00 

46.23 

-13.88 

334.36 

49.41 

-13.88 

5.78xl0-4 

327.96 

41.57 

-10,48 

4 . 82xl0~4 

321.31 

43.43 

-10.48 

4 . 70xl0~  3 

314.93 

53.83 

-13.40 

4.32xl0~3 

308.81 

53.71 

-13.83 

302.94 

54.57 

'-14.18 

297.32 

54.54 

-14.18 

291.94 

53.07 

-14.18 

1.09X10-3 

286.83 

45.35 

-12.58 

281.98 

48.43 

-13,74 

277.41 

50.72 

-15.20 

273.11 

31.93 

-9.95 

269.08 

Tables  3 and  4 present  the  similar  information  for  sensor  5,  We  have 


chosen  to  provide  in  this  report  results  for  one  sensor  from  each  probe 


i 

nn 

r 


assembly  (front  & back)  only.  Calculations,  however,  were  made  for  each  of 
the  sensors. 

In  sunmary:  We  have  a relatively  meaningful  sample  of: 


]±  ..  V ("ake> 

'+  l+  (fr0nt) 


I (wake) 

I+  (ambient) 


which  lends  itself  to  a physical  parametric  analysis  such  as:  a * f(R_,$u) 


Figure  1 shows  the  variation  of 


I (wake) 
a I+  (front) 


with  altitude  using  measure- 


ments from  sensors  1,2,3  (i.e.  from  the  forward  probe  assembly).  It  is  seen 
that  the  variation  in  a extends  over  four  orders  of  magnitude  in  the  altitude 
range  300  km  to  1100  km.  The  steepest  variation  being  in  the  altitude  range 
300  km  to  500  km.  Figure  2 shows  the  variation  of  Rp  ■ f(H).  It  is  seen 
that  the  largest  variation  in  Rjj(-Rq/X  ) is  ln  the  altitude  range  300  km  to 
(450-500)  km.  Figure  3 is  obtained  basically  from  Fig.  1 and  Fig.  2 and 
shows  the  variation  of  a « f(R^).  The  result  obtained  shows  that  the  case 
of  a 'large-body'  becomes  applicable  for  Rp  « 40.  This  is  a new  experimental 
finding  which  is  relevant  to  spacelab  applications  and  is  of  Importance  to 
theoretical  wake  modeling.  In  other  words:  Figures  1-3  show  that  the 
amount  of  ion  depletion  in  the  wake  of  an  ionospheric  satellite  for  ^ 40 
is  several  orders  of  magnitude  larger  than  for  the  cases  were  24  < R^  < 40. 
The  behaviour  of  • f(H)  in  the  altitude  range  600  km  to  700  km  requires 
further  study. 


■ J • 


We  have  repeated  the  same  procedure  for  the  measurements  made  by  sensors 


5,  6,  and  8.  As  seen  from  Figure  4 (u=f(H))  the  data  is  of  lower  quality 
compared  with  the  data  from  sensors  1,  2,  and  3 but  still  depicts  qualitatively 
the  general  behaviour  as  seen  in  Figure  1. 

We  have  also  plotted  [M+)AV=f(H)  (Figure  5)  in  order  to  study  a = f(SAV) 
where:  SAy  = average  ionic  Mach  number.  No  clear  correlation  can  at  this 

point  of  the  study  be  pointed  out.  We  have  also  examined  a = f($N>.  As  seen 
from  Figures  6 and  7 no  obvious  correlation  between  a and  is  obvious.  In 


fact  from  Figure  6 we  see  that  valjes  of  « which  differ  by  about  two  orders 
of  magnitude  have  similar  values  of  4>N  (i . e . , 12 . 5<  |$N|<  15)  which  can  point 
to  the  fact  that  another  parametric  dependence  is  dominant.  And  indeed 
based  on  the  present  sample  of  measurements  we  believe  that  a = f(Rp)  is 
dominant. 

Note:  The  above  results  (as  an  experimental  finding)  is  publishable. 
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3.  Discussion  of  Results 


A.  About  a: 


Optimally,  one  would  seek: 


I+  (141°<0<144°)‘ 
I+  (41°<0£44°) 


p+  (®=i8o®>i 

“ = [l+  (e=0°)  J 


rather  than:  a = 


as  obtained  from  the  present  preliminary  study.  This 


is  evidently  a limitation  of  the  present  data.  Other  data  which  are  now 

being  examined  show  availability  of  measurements  for:  0 . -*•  0°  and  6 

min  max 

180°.  In  any  case  it  should  be  possible  to  apply  the  appropriate  sin/cos 
function  to  I+  (9min)  and  obtain  I+(6=0°)  but  it  is  not  possible  to  obtain 
I+  (0  = 180°)  from  I (0  'o  144°).  This  is  so  since  at  the  present 
time  there  is  no  reliable  theoretical  model  for  the  wake  zone  which 
could  unambiguously  correlate  the  amount  of  depletion  at  0^  with 
the  amount  of  depletion  at  0^  for  i//k,  k > i,  or  correlate  unambigu- 
ously ion  density  ratios  (for  130° £0<18O°)  with  ion  current  ratios  for  that 
angular  range.  It  is  possible  to  use,  say,  a ’neutral  approximation’  (i.e. 
the  approximation  where  ions  are  treated  as  neutral  even  in  the  very  near 
vicinity  to  the  satellite)  and  extrapolate  I to  the  maximum  wake  position 

fi.  (e.)l  [>l+‘Vl 

or  correlate  1^ — jr-  H with  L^q  )j  f°r  and  02  in  the  range  ^O0  180 


but  then  it  would  be  physically  incorrect  to  use  such  results  in  a theory- 
experiment  comparison. 


B.  About  a 


The  sample  of  data  used  (Orbit  1131,  date:  02/24/76)  has  demonstrated 
the  dependence  of: 

fl . (wake)  "1  [i . (wake)  "1 

f<V 


I (wake) 

+ X 

I (wake) 

I (front) 

I+  (ambient) 

In  the  ranges: 

24<R])<56 

10<|^18 

under  the  limitation  that  0 = // 0 and  0 = //180.  Namely,  the  amount  of  ion 

min  max 

current  depletion  in  the  wake  does  not  actually  depict  the  situation  for  the 
maximum  rarefaction  zone.  Nevertheless,  the  study  already  performed  yielded 
new  information  not  yet  published.  It  appears  from  the  present  study  that  a » 
f(Rp)  is  more  significant  than  a possible  a=f(<J>N)  dependence.  It  should, 
however,  be  noted  that  such  a result  does  not  necessarily  have  to  hold  for 
plasmas  specified  by  other  ranges  of  plasma  and  body  parameters.  The  continu- 
ation of  the  study  is  expected  to  contribute  to  the  latter,  and  further  explore 
a = f(SAV)  for  other  flow  regimes. 

The  results  already  obtained  will  have  now  to  be  further  ’refined'  and 
compared  with  other  complementary  results  from  in-situ  and  laboratory  simula- 
tion studies  as  well  as  compared  with  some  theoretical  wake  models.  There  is 
doubt  whether  the  level  of  funding  available  to  us  will  suffice  to  perform 
the  latter.  While  there  is  doubt  whether  a meaningful  theory-experiment 
study  can  be  performed  within  the  alloted  budget,  we  should  be  able  to  con- 
tinue the  experimental  investigation  (via  data-analysis)  and  extend  it  to 

* 

be  applicable  and  relevant  to  future  spacelab  missions. 
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